The Optical Depth to Reionization as a Probe of Cosmological and 

Astrophysical Parameters 



Aparna Venkatesan 

Department of Astronomy and Astrophysics, 
5640 S. Ellis Ave., University of Chicago, 
Chicago, IL 60637, USA 

O 

£S) ' aparna@oddjob.uchicago.edu 

a" 

ABSTRACT 



o 



ON 

Or 



Current data of high-redshift absorption-line systems imply that hydrogen reioniza- 
tion occurred before redshifts of about 5. Previous works on reionization by the first 



J> ' stars or quasars have shown that such scenarios are described by a large number of cos- 



mological and astrophysical parameters. Here, we adopt a semi- analytic model of stellar 



reionization in order to quantify how the optical depth to reionization depends on such 



parameters, and combine this with constraints from the cosmic microwave background 



ON , (CMB). We find this approach to be particularly useful in alleviating the well-known 



degeneracy in CMB parameter extraction between the optical depth to reionization and 
the amplitude of the primordial power spectrum, due to the complementary information 
from the reionization model. We also examine translating independent limits on astro- 
i—i ■ physical parameters into those on cosmological parameters, or conversely, how improved 

determinations of cosmological parameters will constrain astrophysical unknowns. 
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s_j ■ 1. Introduction 

a, 

Observations of the spectra of distant quasars and galaxies have revealed the absence of a Gunn- 
Peterson trough, implying that the intergalactic medium (IGM) was highly ionized by redshifts of 
about 5. Since the universe recombined at redshifts, z ~ 1100, the IGM is expected to remain 
neutral until it is reionized through the activity of the first luminous sources. At present, it 
appears that the reionization of hydrogen occurred before z ~ 5 (Schneider et al. 1991; Hu et al. 
1999; Spinrad et al. 1998), while that of doubly ionized helium is thought to have occurred before 
z ~ 3 (Hogan et al. 1997; Reimers et al. 1997). 

The ionizing sources responsible for reionization can be a variety of astrophysical objects, and 
much work has been done on reionization by the first stars (Haiman & Loeb 1997; Fukugita & 
Kawasaki 1994), the first quasars (Haiman & Loeb 1998a; Valageas & Silk 1999), proto-galaxies 
(Cen & Ostriker 1993; Gnedin 2000; Giroux & Shapiro 1996; Ciardi et al. 2000; Miralda-Escude 
et al. 2000; Madau et al. 1999), and related phenomena such as supernova-driven winds (Tegmark 
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et al. 1993), and cosmic rays (Nath & Biermann 1993). Quasars are natural candidates for both H 
I and He II reionization, as they are more luminous and provide harder ionizing radiation relative 
to stars, and they are seen up to the highest redshifts of current observations. However, there are 
indications of a turnover in the space density of the QSO population, which apparently decreases 
beyond z ~ 3. As this is based on optical surveys, this observed decline is subject to the effects 
of any dust obscuration along the line of sight; however, recent radio observations also appear to 
indicate a declining QSO population beyond z ~ 3, so that the comoving emission rate of Lyman 
continuum photons from QSOs is deficient by a factor of ~ 4 relative to that required for reionization 
(see, e.g., Madau (1998) and references therein). If this is real, then QSOs may be less plausible 
as sources for H I reionization (see also Rauch et al. (1997)). The alternative is to allow QSO 
formation in collapsing halos from the outset, and to postulate a large population of faint QSOs at 
z > 5, with the observed turnover being true for bright QSOs only. 

Stellar reionization is attractive for several reasons. The first stars are expected to form at 
z > 10, and are capable of ionizing hydrogen. Furthermore, they create heavy elements, and may 
account for the low (0.01 Z & ) but persistent metallicity seen in the Ly-a forest clouds out to z ~ 4. 
Finally, the detections of a He II absorption trough in high-z quasar spectra appear to indicate 
a soft component to the UV background (Haardt & Madau 1996), consistent with the ionizing 
spectrum of stars. 

Early work on hydrogen reionization (Arons & Wingert 1972) described the appearance of 
the first luminous sources, about which ionized bubbles gradually expand into a neutral IGM; 
eventually such H II regions overlap and the universe becomes transparent to Lyman continuum 
photons. In principle, only one ionizing photon per neutral atom in the IGM is required, but 
the effects of recombinations ensure that for an atom to stay ionized, the rate of ionizing photons 
generated by sources must be greater than the rate of recombination at that epoch. This is of 
particular importance at high z's, when the IGM had greater density. Just how much more than 1 
photon per baryon is needed is a function of the cosmology and the assumptions of the reionization 
treatment; some evolution with redshift is also expected. A qualitative assessment may be made 
however; for example, from arguments of producing the average IGM metallicity in C IV of about 
10~ 2 Z© at z < 5, Miralda-Escude & Rees (1997) arrived at a requirement of 20 ionizing photons 
per baryon. For the stellar reionization model that will be considered in this work, we will see below 
that about 5 ionizing photons per baryon are available for H I reionization and prove sufficient. 
Not more than a few percent of all baryons need participate in early star formation for reionization 
to occur by z ~ 5, though this number may reach values of up to about 15% (§2). 

Though reionization by early stars would occur at redshifts well beyond current observations, 
it has many distinct consequences which can be feasibly constrained by current and future exper- 
iments. Some of these include, as mentioned above, the evolving IGM metallicity and the cosmic 
ionizing background as derived from spectral features in high-z absorption-line systems. High-z 
reionization will also leave a signature in the CMB through the Thomson scattering of CMB photons 
from free electrons (Sugiyama et al. 1993; Dodelson h Jubas 1995; Tegmark et al. 1994; Tegmark 



-3- 



& Silk 1995; Hu 2000; Zaldarriaga 1997). Depending on the epoch and degree of reionization, we 
expect an overall (somewhat scale-dependent) damping of primary temperature anisotropies in the 
CMB, the generation of new temperature anisotropies on the appropriate scales through the effects 
of second-order processes and the degree of inhomogeneity in the reionization process (Gruzinov & 
Hu 1998; Knox et al. 1998), and finally, the creation of a new polarization signal, as the process 
of Thomson scattering introduces some degree of polarization even for incident radiation that is 
unpolarized. Scattering from the ionized IGM, or the reprocessing of starlight into far-infrared 
wavelengths by dust formed from early supernovae (SNe), will also cause the CMB to undergo 
some spectral distortion (Loeb & Haiman 1997); this can be measured experimentally through the 
Compton y-parameter. These and other observational signatures that have the potential to con- 
strain the epoch, and hence possibly the source, of reionization have been examined in the literature 
(see Haiman & Knox (1999) for a summary). 

A model of reionization is therefore, in principle, eminently testable. Current detections of 
the first Doppler peak in the CMB's temperature anisotropies limit the total optical depth to 
electron scattering, r e , such as may arise from reionization, to be r e < 1 (Scott et al. 1995). 
Future experiments such as the Next Generation Space Telescope (NGST) or the Space Infrared 
Telescope Facility (SI RTF) may detect the high-z sites of reionizing sources (see, e.g., Haiman &: 
Loeb (1998b)), or at least exclude currently viable candidates, while upcoming CMB experiments 
such as MAP or the PLANCK surveyor can measure r e to very high accuracies by combining 
information from temperature anisotropies and polarization in the CMB. 

The optimistic prospects for testing reionization and the increasing multiwavelength view of 
the high-z universe have generated a large body of work on reionization models in the last few 
years, whose techniques fall broadly into numerical (Gnedin 2000; Chiu & Ostriker 2000; Gnedin 
& Ostriker 1997; Cen & Ostriker 1993) or semi-analytic methods (Haiman & Loeb 1997, 1998a; 
Tegmark et al. 1994; Valageas &: Silk 1999). The former have the advantage of being able to track 
the details of radiative transfer, incorporating the dumpiness of the IGM and the essentially non- 
uniform development of ionizing sources, and, perhaps most importantly, describing the process of 
reionization in a quantitative fashion. The advantage of semi-analytic approaches is their inherent 
flexibility and ability to probe the parameter space of a reionization model at will, which is of value 
given the many input cosmological and astrophysical parameters involved. 

For astrophysical sources, the process of reionization is strongly related to the evolution of 
structure in the universe, and could result in feedback effects for subsequent object formation (see, 
e.g., Ciardi et al. (2000)). Of the current theories of structure formation, variants of the standard 
cold dark matter (sCDM) model are considered to be relatively successful at describing the observed 
universe. This picture postulates a critical density universe, with cold dark matter dominating 
the matter content; structures, made up of baryons and CDM, originated in primordial adiabatic 
fluctuations and evolved subsequently through gravitational instability. Current modifications to 
this paradigm include, e.g., the addition of a cosmological constant. 
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The sCDM model, and its variants, are described by a set of parameters which characterize the 
primordial power spectrum of fluctuations, the cosmology of the universe, and quantities related 
to primordial nucleosynthesis. At present, the extraction of such parameters from observations has 
been made feasible by the quality of data from large-scale structure surveys, from cosmic velocity 
flows (Zehavi & Dekel 1999), from Type la SNe (Tegmark 1999), and from current and projected 
future data from the CMB (Zaldarriaga et al. 1997; Eisenstein et al. 1999). Typically, a 9-13 
parameter set describes the adiabatic CDM model, and can be solved for given the data. One of 
these parameters is r e , which is by nature somewhat unique, in that it is the only quantity that is 
not set purely by the physics prior to the first few minutes after the Big Bang. Thus it can po- 
tentially provide clues on post-recombination astrophysics, assuming that the other (cosmological) 
parameters which also affect r e are comparatively well-constrained. 

Several of the semi-analytic works on reionization mentioned above have explored the effects of 
varying model parameters on r e . Other authors have performed CMB analyses that have revealed 
inherent degeneracies in constraining specific combinations of parameters, e.g., r e and the amplitude 
of the primordial power spectrum, A (see, e.g., Zaldarriaga et al. (1997); Eisenstein et al. (1999)). In 
this paper, we examine the results of cross-constraining the range in a cosmological parameter, given 
the allowed band in r e from a reionization model due to astrophysical parameter uncertainty, with 
the permitted range from CMB observations. Specifically, we find that for the combination r e — A, 
the well-known degeneracy in their effects on the CMB can be broken when used in combination with 
the constraints from a reionization model. As r e depends on both cosmological and astrophysical 
parameters, however, such an analysis can be extended to mutual constraints involving these two 
independent classes of parameters, by eliminating r e . The advantage of this is that, given a model 
of structure formation and a reasonable framework describing reionization, as well as the data 
from the CMB, we can use known astrophysics to further constrain cosmology and place tighter 
limits on cosmological parameters, even those that will be determined to high accuracies by future 
experiments. Conversely, a well-determined cosmological parameter can be used to constrain the 
astrophysics of ill-known details of early star formation. This will, at the least, be a powerful test 
of the cosmology, if our understanding of reionization is reasonably correct; the additional hope is 
that this will prove to be an alternative way of constraining the activity of the first stars. 

The plan of the paper is as follows. In §2, we outline the stellar reionization scenario that we 
consider, and set up a parameter set that describes reionization. In §3, we review the standard 
methodology related to parameter extraction from the CMB, and incorporate the parameter set 
from §2 into this formalism. In §4, we present our results, and show the combined constraints from 
a reionization model and the CMB. We present our conclusions in §5. 
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2. The Stellar Reionization Model 



We assume a sCDM primordial power spectrum of density fluctuations, given by P{k) = 
Ak n T 2 (k), where A and n are, respectively, the amplitude and index of the power spectrum, and 
the matter transfer function T(k) is taken to be the form given in Bardeen et al. (1986), with the 
baryon correction as given by Peacock & Dodds (1994). Here, we evaluate the power spectrum 
normalization A through the value of the rms density contrast over spheres of radius 8 h^ 1 Mpc 
today, eg. The cosmology of our model is also described by Qq, the cosmological density parameter, 
Qb, the density parameter of baryons, and h, the Hubble constant in units of 100 km s _1 Mpc -1 . 
We will assume that Oo = 1 throughout this work, and thus will not include it in the parameter 
set to be varied in what follows. We assume that there are no tensor contributions to P{k), and 
set the cosmological constant to be zero. Thus, our cosmological parameter set is [A, n, fij,, h]. 

The reionization of the universe by the first generations of stars is described by the model 
developed in Haiman & Loeb (1997) (henceforth HL97), with the minor modifications described 
below. Briefly, the fraction of baryons in collapsed dark matter halos, Fb, is followed using the 
Press-Schechter formulation; of these baryons, a fraction /* cool and form stars in a Scalo initial 
mass function (IMF) . A fraction f esc of the generated ionizing photons is assumed to escape from 
the host object and propagate isotropically into the IGM. One can then solve for the size of the 
ionized regions associated with each such star-forming cloud, which, when integrated over all haloes, 
yields at each z the average ionization fraction of the universe, given by the filling factor of ionized 
hydrogen by volume (Fjjii)- We assume that the IGM is homogeneous, in which case the ionized 
region created by each source can be taken to be spheres of radius r\. Reionization is defined to 
occur when Fjjii = 1. The total optical depth for electron scattering, T re i on , to the reionization 
redshift z re i on , is given by integrating the product of the electron density, the ionization fraction, 
and the Thomson cross section along the line-of-sight from the present to z re ion- The cosmology of 
the universe enters T re i on through the first two terms of the integrand, and also through the path 
length of the photons last scattered at z re i on . 

Our adopted model is summarized by the following equations for an S7q = 1 universe: 



Fhii(z) = Pb (z) dz on -j^(z on ) (—rf(z on ,z)) (2) 

p Zreion 

T reion = 0.053 n b h dz VTTz' (1 - f*F B (z))F HH (z) (3) 

J 

The critical overdensity S c = 1.686, a(R,z) is evaluated with a spherical top-hat window 
function over a scale R oc Mc, where Mq is the minimum halo mass that collapses at a given 
redshift. While a natural choice for Mq is the baryonic Jeans mass, given by ~ 1O 6 M [(1 + 
z)/100] 1 - 5 , this assumes that collapsing halos at a given mass scale are equivalent to star-forming 
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halos. However, several authors (Tegmark et al. 1997; Haiman et al. 1996) have argued for a higher 
value of M<7, based on the requirement of an effective coolant for baryons in a halo to fragment 
into stars. The picture is as follows: the very first stars form from metal-free gas and cool through 
primordial molecular hydrogen. The universe at that epoch is transparent to photons in the energy 
range 11.2-13.6 eV, but is opaque to more energetic photons. This initial trace level of stellar 
activity easily photodissociates the remaining H2 in the universe (whose abundance relative to H I 
is very low) , well before the H I ionizing flux has built to values sufficient for reionization (see HL97 
and references therein). Subsequent halo formation continues but star formation halts as there is 
no coolant available, and can only resume when halos more massive than 1O 8 M0[(1 + z)/10] -1 " 5 
collapse, which can utilize line cooling by atomic H. We set Mq to be this higher value, with the 
understanding that Fb now represents the fraction of all baryons that are in star-forming halos. 
Finally, z± is the earliest redshift at which the first stars can form, and here we set it to be 100. 
Prior to this, the temperature of the IGM is coupled to that of the CMB, and the CMB photons 
are still energetic enough to photodestroy H~ , thus preventing the formation of H2 through the H~ 
channel, an important cooling mechanism for structures at these high redshifts to fragment into 
stars. 

The evolution of an individual ionization front is characterized by the ionization radius rj, 
and for a time-dependent source luminosity, can be solved for through a differential equation as 
in HL97, where the rate of emission of ionizing photons from a stellar population of metallicity 
Z = 10~ 4 Z Q was constant for about 2 million years before declining with the death of the massive 
stars in the IMF. In this work, we use the analytic solution from Shapiro & Giroux (1987) (SG87 
hereafter) for the evolution of rj in an expanding universe in units of the Stromgren radius r s . 
The Stromgren radius represents the equilibrium reached in the IGM between a source's ionizing 
photon rate and the IGM's recombination rate; r s increases with decreasing z, or equivalently, with 
decreasing average IGM density. The maximum value that T{ can have is r s ; only sources at very 
high redshifts (~ 100) have ionized regions that fill their Stromgren spheres [SG87]. Note that as 
the SG87 solution does not account for time-varying sources, we expect T re i on to be overestimated 
(reionization occurs earlier) compared to HL97, but as we will see in the next section, this is 
a very slight effect. Thus, in this work, vi = (ri/r s )sG87 r s, where r 3 = 3S(0)/(47raBn H (z)), 
«b = 2.6x 10 -13 cm 3 s -1 , and the initial emission rate of ionizing photons leaving the host object 
is 5(0). Fhii{z) (see eq. [2]) is determined at each redshift z by integrating over the product of the 
rate of new halos that formed stars at a turn-on redshift z on (where z < z on < z*), and the ionized 
volume per unit mass associated with such objects, for a source mass M. A detailed treatment of 
the evolution of Fb and Fjjii with z, given various choices of input parameters, may be found in 
HL97. For their standard model, Fb rises rapidly from values of ~ 10~ 3 at z ~ 20 to about 0.1 at 
z ~ 10; during this period, Fun rises steeply from about 10 -4 [z ~ 32) to unity at z ~ 18, so that 
reionization occurs relatively quickly with the growth of Fb- 

We now see that r re j on is a function of several cosmological and astrophysical parameters. 
The astrophysical parameter set is [/*, Mc, 5(0)]. 5(0) is itself a function of several variables: the 
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choice of the IMF, the metallicity Z of the progenitor stars, f esc , /*, and the halo's mass. The last 
two factors account for the fraction of star-forming baryons in each halo, while f esc represents the 
loss of ionizing photons to the host cloud before reaching the IGM. Let us now consider the first 
two variables. Though there have been some arguments for the IMF to be biased towards high- 
mass stars in the early universe (Larson 1998), the details of the nature of star formation under 
those conditions are still not well understood. In the absence of a convincing theory of primordial 
star formation, the most reasonable assumption is to take a present-day IMF and calculate the 
luminosity expected from metal-poor stars. As reionization is affected primarily by the massive 
stars in any IMF, an IMF in the past biased towards high-mass stars would still have the same 
emission spectrum of ionizing photons, while one dominated by low-mass stars is unlikely to reionize 
the universe by z ~ 5. We therefore take S(0) to be as given in HL97 for Z = 10~ 4 Z Q stars from 
standard stellar evolutionary models; it includes the ionizing radiation from stars only, and is steady 
at / + / esc 10 46 photons s _1 Mq 1 for about 2 million years before declining rapidly, which is consistent 
with the value in, e.g., Ciardi et al. (2000). As a rough estimate, this translates to ~ 5 ionizing 
photons per baryon in the universe, for f+ = 0.05, f esc = 0.2. 

The ionizing photon contribution from SNe is relatively small (HL97), but most of the mass 
assigned to forming stars is eventually returned to the IGM. Therefore, the second factor on the 
RHS of equation (3) for T re i on should have an extra contribution /snI^Fe to account for the extra 
baryons that are available for new star formation, where fsN is an IMF-averaged fraction of the 
progenitor mass that is expelled into the IGM at the end of the star's life. We expect between ~ 
50-95% of the parent star's mass as ejecta, until a point is reached (for stellar masses ranging from 
50-100 Mq) where the entire star collapses into a black hole. For the low values of that will 
be considered here, the product f±fsN will be a small correction and may be ignored. Moreover, 
the mass of the ejecta from a dying star depends sensitively on the stellar metallicity, with low-Z 
stars having higher remnant masses and less ejected material relative to solar-Z stars (Woosley 
& Weaver 1995). Thus, /sn is likely to be highly variable, both spatially and with z due to the 
evolving metallicity of subsequent generations of stars, and is more appropriately modelled in a 
simulation rather than in a semi-analytic model. The calculated values of r re j on here may be taken 
as a lower limit. 

This leaves the astrophysical parameters, f esc , and Mc- We note that in most semi-analytic 
models, /* is set by the choice of Mc, as the stellar Z-output, particularly in 12 C, is combined with 
the evolution of Fb to produce the observationally detected average carbon abundance of 0.01 Z & 
in the Ly-a forest clouds at z ~ 3 (Songaila & Cowie 1996). Thus, /* and Mc are not independent 
of each other if we choose the above normalization; for the HL97 choice of Mc, /* = 0.13. This is 
the maximum value that can have from arguments of avoiding IGM over-enrichment; given the 
approximately order-of-magnitude scatter in the average metal abundance of the Lyman-a clouds, 
and that one need not require the reionizing stars to solely account for Zjcm, /* may be smaller 
than ~ 0.15. 



-8- 



As an aside for the interested reader, we note here two drawbacks of normalizing /* via C. 
One is that the massive stars in the IMF (> 10M Q ) are the ones relevant for reionization, while 12 C 
is produced dominantly by the intermediate-mass stars (2-8 Mq). Thus, if the IMF was different 
in the past, the carbon abundance in the Ly-a forest does not constrain the massive or reionizing 
stellar activity in early halos. The second point to note is that the pause in star formation caused 
by the initial dissociation of H2 led to the choice of Mq, as in HL97. This minimum halo mass 
corresponds to objects of virial temperature of ~ 10 4 K. At this temperature, the host object is 
immune to photoionization heating (as pointed out in HL97), and so if outflows are desired to expel 
the generated metals into the IGM (i.e, the Ly-a forest), the mechanical energy of SNe must be 
invoked. Again, the massive reionizing stars end their lives as Type II SNe an order of magnitude 
in time before their intermediate-mass carbon-producing compatriots. As all the mechanical input 
lies with Type II SNe, the question arises of how the carbon, produced significantly later, leaves 
the host halo to mix with the IGM. Furthermore, Type II SNe occur on much more predictable 
timescales, i.e. immediately following the progenitor's death, than do Type la SNe (3-10 Gyr), 
and there is not much more than a wheeze to be had from the deaths of intermediate-mass stars 
as planetary nebulae. 

Having voiced these objections, we point out that while the 12 C connection as made above 
between Mq and /* is not ideal, postulating a general relation between these two variables is not 
ad hoc. The value of Mq does intrinsically determine the stellar history, metallicity and luminosity 
evolution of the universe; the high value of Mq in HL97 and other works is physically well-motivated 
by the necessary step of having an available coolant to aid star formation. We proceed to set Mq 
= 10 8 M Q [(1 + z)/10]~ L5 for the semi-analytic treatment here, and now narrow our astrophysical 
parameter set to [/ esc , /*]. 

We end here by addressing some of the issues that are not accounted for in this work. The 
IGM is assumed to be homogeneous, but clearly some dumpiness will develop in the IGM from 
the growth of initial density inhomogeneities, and the assumption of the average ionized fraction 
at a given redshift being equal to the H II filling factor will eventually break down. However, 
this appears to be a relevant effect only at "late" times (z < 10), when the fraction of baryons 
in collapsed structures becomes significant (Gnedin & Ostriker 1997), or for baryon-dominated 
universes [SG87]. Therefore, we will assume that the clumping factor is unity (homogeneous IGM) 
for the rest of this work. We have also neglected corrections from doubly ionized helium, which 
is not problematic as the spectrum of photons produced by stars is softer than that from quasars, 
and is more relevant for H I rather than He II reionization (see, however, Tumlinson & Shull (2000) 
on the helium-ionizing spectrum from zero-metallicity stars). We have set Fh £ jj = Fhjj, but this 
introduces an error of not more than a few percent (Tegmark & Silk 1995). Finally, we have not 
included the effects of bias in the normalization of the matter power spectrum, i.e., we assume that 
light traces the underlying mass distribution. 
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3. Constraints From The Microwave Background 



As discussed in the introduction, signatures from reionization are expected in the CMB; 
an accurate measurement of T re i on or the detection of post-recombination features in the CMB 
anisotropies have the power to constrain the reionization epoch and the nature of the sources 
through the angular scale #(oc /~ 1 ) on which they affected the CMB. Here, I is defined from ex- 
panding the angular power spectrum of the CMB in terms of its multipole moments Ci and Legendre 
polynomials: 



The effect of r re j ori is to introduce an overall damping of the temperature Cis (HL97, and 
references therein), except at the largest scales. As discussed by, e.g., Zaldarriaga et al. (1997), 
this is practically indistinguishable from the generally reduced values of C/ expected from simply 
having a lower amplitude, A, of the primordial power spectrum; the difference between these two 
effects at the smallest /'s is obscured by cosmic variance. While the amount of the damping due 
to 

Treion is /-dependent and can potentially break this degeneracy, the accuracy to which T re % on can 
be estimated from temperature anisotropy maps alone is not sufficient to distinguish between these 
two different effects (Zaldarriaga 1997). 

When combined with the polarization data from the CMB however, r re j on can be constrained 
with far greater precision. Linear polarization is generated by the primary temperature quadrupole 
anisotropy photons scattering off the free electrons in the reionized IGM, and is a relatively clean 
probe of the epoch of reionization (Zaldarriaga 1997). The polarization signal is expected at 
low levels compared to that from temperature anisotropies, and may prove difficult to measure, 
especially for low optical depths. Nevertheless, r re j on can be detected, in principle, by future 
experiments to within 1-sigma errors of, e.g., 0.69 (0.022) without (with) polarization information 
for MAP, and 0.59 (0.004) correspondingly for PLANCK (Eisenstein et al. 1999). 

We wish to combine the constraints on cosmological or astrophysical parameters from a reion- 
ization scenario with those from the CMB; in order to do this for the latter, we follow the standard 
prescription as outlined in, e.g., Jungman et al. (1996) and Knox (1995). We assume Gaussian 
initial perturbations, and that the multipole moments C\ are determined by a "true" set of N 
theoretical parameters [P/v]. If we define the likelihood function C of observing a set of C/S, given 
P/v , then the behavior of C near its maximum can be quantified in terms of the Fisher information 
matrix, whose elements are given by the second derivative of the logarithm of C with respect to 
pairs of parameters in P/v- The Fisher matrix then represents the accuracy with which P/v can be 
estimated from a given data set, here the CMB's experimentally measured C/s. Further assuming 
that C has a Gaussian form near its maximum, the Fisher matrix is given by: 




(4) 



1=2 
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(5) 



where ci; is a measure of how the observed C;s are distributed about the mean value of the 
true Ci(Pn)s. We assume that o\ is cosmic- variance-limited, and ignore terms arising from the 
instrumental noise associated with an experiment and from any foregrounds. For a sky fraction 
/sky that has been mapped, o\ can be approximated by: 



) 



° l ~ \J (2/+l)/ Bky °l( P N) ( l < l max 
OO (I > l ma x) 



(6) 



We will consider two cases here: (/ ma;r . = 400, / s k y = 0.01), which is roughly representative of 
data from current CMB experiments, and (l ma x = 3000, / s k y = 0.8) for the data expected from 
PLANCK . As we neglect any experimental or systematic effects, the power of the C/s to constrain 
P/v, as presented here, is the "best possible" case. Note also that the above formulae are valid 
when only the temperature information from the CMB is used. More general expressions for the 
case of including polarization data may be found in, e.g., Zaldarriaga et al. (1997). 

The derivatives of the Qs with respect to Pjy were computed for each parameter using two-sided 
derivatives with step sizes being chosen so that the value of this derivative remained stable (see, e.g., 
Eisenstein et al. (1999), Appendix B.l). The C;s themselves for a given parameter set were found 
using the publicly available CMBFAST (v. 2.4.1) 1 . Note that the parameter set describing the 
reionization model is [A, h, n, f esc , /*], which yields r re j ori , whereas the CMB data can determine 
the cosmological parameters and T reion , or equivalently, [P C osmo, r re ion(Pcosmo, -Pastro)]- Therefore, 
when we specify r (e.g, to CMBFAST), the cosmological and the astrophysical parameter sets 
(-Pcosmoj Pastro) are no longer independent, but are related through the reionization model, and the 
C\ derivatives become: 



dd 



dP n 



dCi 
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dCi 
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dQ 
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dQ dr 
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(7) 
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Once the Fisher matrix Fij has been constructed, it can be inverted to give the covariance 
matrix C between the parameters P/v; Ca represents the minimum variance in the estimate of Pj. 
Any 2x2 submatrix of C can then be extracted, giving the ellipse equation for the joint confidence 
region in the 2-parameter subspace of interest, 



1 http://www. sns.ias.edu/~matiasz/CMBFAST/cmbfast. html 



[P-P n ].(C 2x2 )- 1 .[P-Pn] = A, 

where A is set throughout this work to be at the 68% confidence level. 



(9) 



4. Results 

We present our results here from combining the reionization model (§2) and the constraints 
from the CMB (§3). This analysis assumes that the density perturbation spectrum at the CMB 
and structure formation scales is described by the same power law. For the choice of parameters 
in HL97, where /* = 0.13 and f esc = f esc (z), 

we obtain T re j on — 0.0734, Or Zreion ^ 

18.4, which 

only slightly exceeds the HL97 value of r re j on = 0.07. Henceforth, we will refer to r re j on as r for 
convenience, noting that r is always evaluated to the reionization epoch in this work. We define 
our standard model (SM), fixing Q = 1, as given by A(a s = 0.7) = 1.55 xlO 6 , Q b = 0.05, h = 0.5, 
n = 1.0, f esc = 0.2, /* = 0.05, and r = 0.0573, with reionization occurring at z ~ 15.5. 

1.0 

0.5 



h 0.0 



-0.5 
-1.0 

1.0x10" 2.0xl0 u 3.0x10" 4.0x10" 5.0x10" 

A 

Fig. 1. — Combined constraints in the t-A plane from the reionization model and current CMB 
data. The standard model (SM) is given by the solid line, the shaded band represents the as- 
trophysical uncertainty in the reionization model, given the allowed ranges for (fescf*), and the 
ellipse is the 1-a joint confidence region from current CMB data. 

As a simple example, we begin with the t-A plane, shown in Figure 1, where we isolate the 
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A 



Fig. 2. — Magnified version of Figure 1; the additional light nested band represents the uncertainty 
in the value of f esc alone. 

dependence of r on A, keeping all the other parameters in the SM fixed. The range of ag is ~ 
0.5-0.8, from the large-scale distribution of clusters of galaxies (see, e.g., Bunn &; White (1997), 
and references therein), but is 1.2 when normalized to COBE for the SM choice of cosmological 
parameters. As an illustrative range for the plots, we normalize A for a$ = 0.5-1.2. The solid line 
represents the SM, while the light shaded region represents the uncertainty due to the astrophysics 
of the reionization model (or AS for astrophysical slop), with [f esc , /*] = [1.0 (0.1), 0.15 (0.01)] for 
the upper (lower) envelope. The range for /* (0.01-0.15) is set as follows: the lower limit comes 
from numerical simulations of star formation (see, e.g., Ciardi et al. (2000), and references therein), 
while the upper limit corresponds approximately to that in HL97 and Haiman & Loeb (1998a). The 
value of f esc has been estimated through a number of theoretical and observational methods (see 
Wood h Loeb (2000), and references therein). Here, we take the range for f esc to be 0.1-1.0, the 
lower limit coming from Dove et al. (2000) (henceforth DSF), who modelled the escape fraction of 
ionizing photons from OB stellar associations in the H I disk of the Milky Way, and found that for 
a coeval star formation history, f esc = 0.15 ±0.05. Our choice of this limit from DSF is motivated 
by the similarity of their model's luminosity history to that in HL97; as noted above, there are 
alternate values for f esc in the literature for a variety of astrophysical environments. We will use 
both the full range for f esc and the more narrow DSF band in later plots. 

To combine this with the CMB constraint, a shaded ellipse representing the 2-parameter 68% 
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Fig. 3. — Constraint from current CMB data in the f esc -A plane, as extended from Figure 2; the 
horizontal band is the permitted range of 0.1-0.2 for f esc from DSF. 

joint confidence region is overplotted, assuming that the true model describing the universe is given 
by the SM and [r, A] = [0.057, A(a$ = 0.7)]. This ellipse is narrow enough that we show a magnified 
version of Figure 1 in Figure 2, with an additional lighter band, nested within the AS band, showing 
the effect of varying only f esc while fixing f+ to its SM value. We see that the CMB does not really 
constrain r or A separately at all, a near-degeneracy that was expected from the discussion in the 
previous section. However, the combination of the CMB confidence region and the AS band is 
much more constraining: this translates to a 1-cr error of about 0.02 for as, which is noticeably 
better than the corresponding value of ~ 0.2 from the CMB ellipse alone. 

We now extend Figure 2 to connect, through r, two a priori independent parameters, f esc 
and A, shown in Figure 3. The purpose of this plot is to probe the potential of cosmology and 
the astrophysics of reionization to constrain each other, given a reionization model. The ellipse 
here, as in Figure 1, reveals the inherent degeneracy between r-related quantities and A through 
the long narrow ellipse. We now overplot the DSF permitted band for f esc (0.1-0.2); clearly even 
this approximate range in f esc considerably narrows the allowed range in A. We note again that 
our choice of the DSF range for f esc was motivated for the reasons outlined earlier, and that other 
ranges for f esc are possible; the main point that is demonstrated by Figure 3 is the power of using 
any such band of independently known astrophysics to constrain a cosmological parameter. Note 
also that A would have to be known to great accuracy to place any limits on f esc that are stronger 
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Table 1: V AND DIAGONAL ELEMENTS OF W from F = UWV T 



4 


-0.3922 


0.1470 


-0.5116 


0.7498 


-0.0277 


0.000 


n b 


-0.5024 


-0.8592 


0.0916 


-0.0320 


-0.0028 


8.42E-16 


h 


0.5012 


-0.3673 


-0.7611 


-0.1856 


-0.0111 


3.613E-15 


n 


-0.5840 


0.3242 


-0.3857 


-0.6342 


-0.0545 


1.599E-14 


fesc 


0.0098 


-0.0039 


0.0111 


0.0041 


-0.2539 


0.9671 


f* 


0.0373 


-0.0148 


0.0420 


0.0155 


-0.9652 


-0.2544 


W 


18110.877 


9251.715 


330.074 


12.455 


0.0191 


6.366E-16 



Note. — Results of the Singular Value Decomposition of Fq X 6 = UWV T ; V, weights w for each column in V, and 
parameters associated with each row in V are shown. 

than the DSF band. 

So far, we have been using the Fisher matrix formalism for specific pairs of parameters, while 
fixing the values of the other parameters in the SM. The more general and proper way to do this is 
to construct a 6 x 6 matrix for the parameter set [A, h, n, f esc , /*], which yields r. However, 
the analysis described in §3 implies that only any five of these six parameters will be independent, 
as the CMB data will determine the cosmological parameters and r. Indeed, the 6x6 matrix, when 
constructed, proves to be singular. We note here briefly some informative aspects of performing 
singular value decomposition (SVD) (Press et al. 1992) of Fq x q, so that F = UWV , where W is 
a diagonal 6x6 matrix containing the singular values w. An element that has an anomalously low 
value (close to zero) in W implies that the corresponding column in V is a linear combination of 
parameters that will not be well-constrained. 

Table 1 shows the matrix V with the corresponding column weights w resulting from this 
decomposition; also shown are the parameters associated with each row in V, with A expressed 
as <7g. We see that the sixth weight is very close to zero, so that the sixth column of V contains 
combinations of [P C osmoj Pastro] that are poorly constrained by this analysis. This column has terms 
corresponding to essentially only Pastro the dominant contribution coming from f esc . This implies 
that f esc will not be well-determined from the CMB (via r), given the reionization model considered 
here. The first five columns of V also convey what combinations of these six parameters will be 
constrained; we note that f esc has very small contributions in these, i.e., to the information that can 
be extracted from the CMB. In comparison, f± can be better determined from the CMB, as seen 
from columns 1-5, particularly the fifth, where the dominant term is from This insensitivity 
of the CMB data to f e sc can be traced back to the stellar reionization model we adopted here; 
variations in /* affect r more significantly than do those in f esc (see Figures 12 and 15 in HL97). 
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The covariance matrix for F can be found by, C = F^ 1 = VW~ X U T . As the ratio of W's 
minimum to maximum value, ~ 3 xl0~ 20 , is very small compared with machine precision, we 
follow the usual technique of adjusting the anomalously low singular value in W, here wq, to zero 
(Press et al. 1992); despite this, the SVD inversion of F still produces an inaccurate covariance 
matrix, i.e., C • F ^ I. 
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Fig. 4. — Constraints from the CMB in the f esc -h plane; larger ellipse represents current data, and 
smaller nested ellipse (line) is from PLANCK. 

We now proceed to work with the independent 5x5 subportions of the full 6x6 matrix, 
which translates to P C osmo and any one of P as tro- These 5x5 matrices are inverted, and the 2x2 
submatrix of interest is projected into the 2-parameter plane as the appropriate error ellipse, which 
displays the confidence region after marginalizing over the other parameters. The results of this 
general Fisher matrix analysis are presented below for the idealized specifications of current data 
and for those expected from PLANCK (§3). Only the temperature anisotropics from the CMB 
are used for Figures 4—11; the polarization information expected from PLANCK is included for 
Figures 12-13. These plots are intended to be examples of the constraints in various Pastro _ -fcosmo 
subspaces. 

Figure 4 shows the case of f esc vs. h; the larger ellipse corresponds to current CMB data, and 
the nested one (appearing as a tiny line) is from PLANCK. For all subsequent cases, we show 
these ellipses separately; the astrophysical range for f esc is omitted from this plot for visual clarity. 
Figure 5 shows the results expected from PLANCK alone for f esc vs. h, with the light shaded 
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Fig. 5. — Constraint from PLANCK in the f esc -h plane; light shaded band represents the entire 
allowed astrophysical range of 0.1-1.0 for f esc , and the dark shaded band represents the permitted 
values of 0.1-0.2 for f esc from DSF. 
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horizontal band representing the full range of f esc (0.1-1.0), and the dark band representing the 
DSF values for f esc (0.1-0.2). The case of f esc vs. is shown in Figures 6 and 7, for current CMB 
data and PLANCK respectively, with the overplotted bands being the same as in Figure 5. 

Figures 8-9 and 10-11 display the respective cases of /* vs. <j| (where erf oc A), and /* vs. 
n. For these 4 plots, the horizontal dark band represents the maximum astrophysical range of 
0.01-0.15 for /*; values below this range are unlikely to be sufficient for reionization and values 
above this range must invoke IMFs other than that of present-day galaxies in order to not violate 
metal production or background light constraints. 

We note here some generic features of Figures 5—11. In all cases, the inclusion of known con- 
straints on the astrophysical parameters strengthens the CMB's limits on cosmological parameters, 
even for the data expected from PLANCK. This is particularly the case with due to r's 
greater sensitivity to /* relative to f esc . Thus the 1-a error for /* from the CMB is significantly 
smaller than that for f esc for all the cases shown here, making independent limits on the former a 
more powerful complementary constraint for cosmological parameters extracted from the CMB. As 
some illustrative examples involving PLANCK data, the entire astrophysical permitted band for 
/* reduces the 1-a error for erg from about 0.02 to less than 0.01 (Figure 9), and for n, from 0.006 
to ~ 0.004. The power to increase such constraints will only become better as /* (or f esc ) become 
better constrained themselves, but it may not matter much for most cosmological parameters in the 
post--PL^4iVC-ftr era (as they will already be determined with great precision), with the exception 
of A or as- For this case, the method here may prove to be a valuable cross-check. 

Figures 5-11 also reveal that even the most promising cases of cosmological parameter determi- 
nation from the CMB's temperature information will not help to constrain astrophysical parameters 
such as /* or f esc , whose currently known ranges as shown through the horizontal bands in each 
figure are typically much smaller than what would be deduced from the joint confidence region. 
This is partly due to the low value of r itself of ~ 0.06 in our SM, which hinders its accurate 
determination from the CMB data. 

When polarization is included for the projected data from PLANCK, we see, from the two 
examples shown in Figures 12-13, that f esc can be determined to about the same accuracy as 
the DSF allowed band, but that the 1-a error for /* is significantly smaller than its astrophys- 
ical uncertainty. Thus, future CMB data may be able to constrain the astrophysical aspects of 
reionization. We recall, however, that we have neglected effects from experimental noise or from 
foregrounds in our analysis, which will enlarge the joint confidence regions in all the figures. While 
this only strengthens the argument of the power of astrophysical limits in constraining cosmology, 
the converse situation, which appears hopeful from Figures 12-13, is realistically tentative at best. 
In short, it is possible in principle that the astrophysics of a stellar reionization model can be 
constrained by limits on cosmological parameters from PLANCK's temperature and polarization 
data, though this may prove difficult to achieve. We may simply have to await the view from 
SIRTF and the NGST to determine the reionizing activity of the first stars! 
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Fig. 7.- 



Same as Figure 6, but for data from PLANCK. 
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Fig. 8. — Constraint from current CMB data in the /*-<t| plane, where erg oc A; shaded band 
represents the permitted astrophysical range of 0.01-0.15 for 



1.0 



0.5 



0.0 




-0.5 



1.0 



0.44 0.46 0.40 0.50 0.52 0.54 

_ 2 



Fig. 9. — Same as Figure 8, but for data from PLANCK. 




Fig. 12. — Constraint from PLANCK in the f esc -^b plane using temperature and polarization; 
lines represent the allowed range for f esc from DSF. 
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Fig. 13. — Constraint from PLANCK in the /*-n plane using temperature and polarization; the 
allowed astrophysical range for /* is the entire y-axis. 
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5. Conclusions 

We have examined the power of a reionization model, given its many cosmological and astro- 
physical parameters, to constrain these input quantities when combined with parameter extraction 
from the CMB. In the case of the well-known degeneracy between r and A in their effects on the 
CMB, we have found that this can be alleviated by the complementary information from a reioniza- 
tion model, and that this remains a useful cross-constraint even when allowing for the astrophysical 
uncertainty in r. 

When we eliminate r and perform a more general Fisher matrix analysis, we find that the 
astrophysical details of reionization can be useful in further constraining the CMB's limits on 
cosmological parameters, even in the case of the expected temperature data from PLANCK . We 
have shown that independent limits on the astrophysical inputs to reionization, despite the current 
uncertainty in their values, reduce the errors for cosmological parameters by a factor of at least ~ 
2. Given that we have considered the most optimistic parameter yield from CMB experiments (§3), 
the use of known astrophysics can only become more valuable for realistic experimental results. 
This is of particular value for as (or A) and n, given their implications for structure formation and 
for theoretical models of the origin of the seeds of structure in the early universe. 

The converse situation- using a projected exquisite determination of a cosmological parameter 
to constrain astrophysical reionization parameters- does not yield quite as interesting results with 
temperature data from current experiments or from PLANCK, even though we made the most 
optimistic assumptions; the 1-a errors for f esc or /* are larger than what are already known to 
be reasonable. When the projected polarization data from PLANCK is included, we found that 
/* in particular may be constrained to far greater accuracy than its current astrophysical uncer- 
tainty; in practice, however, this may prove difficult to achieve, given the effects of foregrounds and 
instrumental noise which we have neglected here. 

In summary, one may take away that the astrophysical details of reionization can strengthen 
the limits on the cosmology of our universe, beyond even the projected parameter yield from future 
CMB data, and that there is more potential to a measurement of r than the determination of 
a single number out of a large parameter space describing adiabatic CDM models. These broad 
conclusions are naturally subject to the assumptions made in this analysis. The sizes of joint 
confidence regions derived from the CMB data for any 2-parameter subspace is determined by 
the full covariance matrix, whose elements' values are dependent on the dimension of the chosen 
parameter space and the selected parameters. The inclusion of more parameters has the generic 
result of increasing the sizes of the error ellipses; therefore, the primary results of this paper can 
only be strengthened when parameter spaces larger than that analyzed here are considered. 

In the sCDM cosmology assumed here, the values of r in our standard model were relatively 
low (~ 0.06). In an open universe, or one dominated by a cosmological constant contribution, we 
expect larger average values of r for a fixed reionization model, as structures freeze out earlier, 
resulting in a longer line-of-sight to the last scattering surface at the reionization epoch. Increased 
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r's can also result from higher values of /* or f esc , or from a lower value of Mc (§2), which would 
allow the first stars to form earlier and more ubiquitously. As higher r's will have a better chance 
of being accurately determined from the CMB, it would be interesting to analyze the constraints in 
this paper, from both the reionization scenario and the CMB, for a more general parameter space; 
we examine this in a forthcoming work. 

This work was presented as part of a dissertation to the Department of Astronomy and Astro- 
physics at The University of Chicago, in partial fulfillment of the requirements for the Ph.D. degree. 
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CMB was performed using CMBFAST by Uros Seljak and Matias Zaldarriaga. This research was 
partly supported by the NSF through the collaborative US-India project No. INT-9605235, by 
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